Orbital ordering and magnetic structures in Sr2-xLaj;FeMo06 and Sr2_xLa2:FeW06 double 

perovskites. 
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We analyzed the possible magnetic and orbital orderings of double perovskites, using a simple 
extension of the double exchange model well suited for these compounds. Orbital ordering is favored 
by the on site repulsion at the Fe ions. We obtain a rich phase diagram, including ferri- and 
antiferromagnetic phases, which can, in turn, be metallic or insulating, depending on the existence 
of orbital order. 
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I. INTRODUCTION 



Ferromagnetic double perovskites of composition 
Sr2_a;La2;FeMo06 (where La can be substituted by other 
double valency ions), known since some time! have at- 
tracted a great deal of attention recently, because of 
their high Curie temperature, metallic character, large 
magnetoresistance and potential applicationsSiSiiiSiSiL^iS. 
While lattice distortions and phonons do not seem to 
play a major role, disorder and in particular the antisite 
defects play a significant role in the properties of these 
compund o^°i-^^i^^i^^i^'^i-^^i^^i^^ . Effects of doping have also 
been extensively studied -^^i-^^i^°i^-^i^^i^'^i^'*i^^ recently . Re- 
lated family of compounds based on Sr2FeW06 are, on 
the other hand, insulating and antiferromagnetioSSiSl. 
Solid solutions including Mo and W ions have also been 

Studied28'29,30^ 

The magnetic and electronic structures of the double 
perovskites with composition Sr2-2:La2,FeMo06 admit a 
simple description*. The Mo^+ and Fe'^+ ions are lo- 
cated at the alternate nodes of a simple cubic lattice. 
The strong exchange coupling within the d orbitals of 
the Fe ion leads to the formation of a spin 5/2 moment, 
from electrons which occupy the exchange-split two eg 
and the three t2g orbitals. The three remaining three t2g 
orbitals at the Fe sites are hybridized, through the inter- 
mediate O ion, with the t2g orbitals of the neighboring 
Mo ions. Because of the symmetry of these orbitals, hop- 
ping between them can only take place along two of the 
three lattice axes, leading to three decoupled two dimen- 
sional bands. Band structure calculations suggest that 
the direct hopping between the Mo t2g orbitals is not 
negligible (this hopping does not change the two dimen- 
sional nature of the bands). The number of electrons in 
the conduction band, per unit cell, is x = 1 -f a; where x 
is the concentration of the divalent La^"*" ions. Most of 
the calculations reported below are given in terms of the 
band filling, x, and the correspondence with the doping 
level, X, is highlighted when needed. 

These arguments suggest that the minimal descrip- 
tion of the electronic structure of these materials requires 
three parameters: the energy difference between the t2g 



levels in the Fe and Mo ions, A = cmo — efc, the Fe-Mo 
hopping, ipe-Mo, and the Mo-Mo hopping, iivio-Mo- We 
assume that the exchange splitting between the spin t2g 
states in the Fe ions is much larger than the width of 
the conduction band. Hence, we need to consider one 
spin state at each Fe ion. This truncation of the states 
used to describe the conduction band is similar to the 
one used in reducing the two band ferromagnetic Kondo 
lattice to the single band double exchange model, used 
in the study of Lai-ajCaa^MnOa and related compounds. 
This tight binding model can also be used to study the 
Sr2FeW06 compound—. 

The phase diagram of the hybridized Mo-Fe system 
described above has been studied using the dynamical 
mean field method^^, and within an exact treatment of 
the electronic wavefunction which allows for the study 
of defects^. Using realistic values for the parameters 
the Curie temperature is comparable to the one exper- 
imentally observed. The generic features of the phase 
diagram have some resemblance to those of double ex- 
change model, to which it can be reduced in some limits. 
One of the most remarkable differences is the appear- 
ance of an antiferromagnetic phase even in the absence 
of direct antiferromagnetic interactionsSS.. 

The analysis outlined above implicitly assumes the 
equivalence of the three t2g orbitals at the Fe and Mo 
sites, and it neglects the possibility of non trivial types 
of orbital order. Electron electron interactions may break 
the symmetry between the t2g orbitals. Such phases 
have been extensively investigated in the related dou- 
ble exchange compounds of composition Lai_2:Caa;Mn03 
and related materiala^^i^'^i^^i^^i'^'''i^^i^^. The electron elec- 
tron interaction in Sr2_a;La2:FeMo06 can be large, as 
the conduction band has a significant weight on the Fe 
orbitals. In the present paper, we analyze the possi- 
bility of the breaking of the orbital symmetry in the 
Sr2-2;La2;FeMo06 family of compounds, using mean field 
theory. The electron-electron interactions are described 
by means of a local repulsion term, U, between elec- 
trons in different orbitals at the same Fe site. Orbital 
ordering has already been reported for the stoichiomet- 
ric Sr2FeW06 compoundSi. We obtain a rich phase di- 
agram, as a function of the number of electrons in the 
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FIG. 1; (Color online). Sketch of the interactions included in 
Hkin, eq.Q. 



conduction band and the difference between the Fe and 
Mo levels. The method is discussed in the next section. 
Then, we present the main results. Section IV discusses 
the relevance of our findings for the physical properties 
of these materials. 



II. THE MODEL AND THE METHOD OF 
CALCULATION. 

We assume that the conduction band arises from the 
hybridization of t2g orbitals at the Mo and Fe sites. 
There are three spin polarized orbitals at the Fe ion, 
and six orbitals at the Mo ion. The hopping takes place 
within planes determined by the symmetry of the orbitals 
(for instance, the hybridization does not allow for hop- 
ping between dxy orbitals along the z axis). There is 
hopping only between states with the same spin. The 
levels at the Fe sites are shifted with respect to those 
at the Mo sites by an energy A. Within these approxi- 
mations, the kinetic energy of the conduction electrons is 
described by three independent hamiltonians of the form: 
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The operators Ca.i,s create electrons with spin s in the 
orbital a at a Mo site, z, and the operators d^ j ^ have 
the same effect at Fe sites. Note that the index s is either 
1 or —1 at the Fe sites, while it can take both values at 
the Mo sites. We have conserved the spin index at the Fe 
sites in order to describe antiferromagnetic phases along 
(111) planes, like that observed in Sr2FeW06 (see below). 
A schematic representation of the couplings in eq.(Q is 
shown in Fig.^. 



The hamiltonian TYkin is defined on a square lattice. 
The resulting energy bands depend on the configuration 
of the core spins at the Fe sites. We will consider the 
two textures shown in Fig.|2]. When all core spins are 
parallel, the unit cell includes one Fe and one Mo site. 
The subbands with spins antiparallel to the core spins 
are built up from one orbital at the Fe site and another 
orbital at the Mo(W) site. The energies, ek, are given 
by: 
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The band with spin parallel to the Fe core spins is local- 
ized at the Mo sites. The energy is: 



Ek = 4tMo-Mo cos{]<ixa) cos(^a) 



(3) 



In the antiferromagnetic configuration shown in the right 
side of Fig. [2] the unit cell includes four sites. The spin 



up and spin down bands are degenerate, although they 
are localized in different regions of the lattice. The spin 
up band is derived from two orbitals at the Mo sites, and 
the orbital at the Fe site whose core has spin down. The 
energies are given by: 
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FIG. 2: (Color online). Configurations of the core spins at 
the Fe ions considered in the text. The rectangles denote the 
corresponding unit cell. 



Det 



iMo-Mo(l + e" 






)(1 



-iU.ya\ 



tMo-Mo(l + e*-")(l + e'^y'^) iFc^Mo(l 



iFc-Mo(l ■ 



iWxa\ 



-fk 

iFc-Mo(l + e 

I — 



■ik^, a \ 



ipc 



-Mo(l 
A„ - 



g-^k^a-) 
iki, a \ 



- e 

Ek 



(4) 



We also include a repulsion between electrons in different 
orbitals at the Fe sites. This term has the form: 
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We use the Hartree-Fock method to analyze the effect of 
this term. Then, assuming that the occupancies of the 
different orbitals at the Fe site are not the same, the Fe 
levels are given by: 
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Here, epc and eMo are defined in the absence of interac- 
tion corrections. The values of the expectation values 
i^ia' s^i,a',s) have to be calculated selfconsistently, in- 
serting the levels given in eq.jnjl into eq.Q. 

Note that the filling of the conduction bands, x can 
vary between and 9, < x < 9. In the following we 
will show calculations for band fillings within this range, 
although fillings such that x > 2 cannot be obtained by 
substituting the divalent Sr ions by trivalent ions. 

From the total energy versus doping curves (see below) 
we identify regions of phase separation in the phase dia- 
gram using the Maxwell construction. This construction 
allows us to identify densities at which different phases 
can coexist as they have the same chemical potential. 



III. RESULTS. 



In the following, we use as unit of energy ipe-Mo- 
From band structure calculations^, this parameter is 
^Fc-Mo ~ 0.3eV. We also fix iMo-Mo = 0.2iFo-Mo and 
U — 12iFo-Mo- These values are less well defined, 
but consistent with existing band structure calculations^. 
The results do not change qualitatively for other values 
of these parameters, although the finer details are modi- 
fied. We vary A = cfc — cmo + U/3 from (strongly hy- 
bridized band, probably adequate for Sr2_j:Laa;FeMo06) 
to A = — 4iFo-Mo (conduction electrons mostly localized 
on Fe ions, valid for Sr2_sLaa;FeW06). The filling of the 
conduction band can vary from to 9. Experimentally 
relevant values of band filling (x) lie between 1 and 2, 
corresponding to doping levels, x, such that < a; < 1. 

Results for the total energy of selfconsistent solutions 
for A — — 4iFc-Mo are shown in Fig.jS]. This case 
favors the localization of the conduction electrons at 
the Fe sites, a situation which probably describes well 
Sr2_2:La2:FeW06^^. The energy difference (per site) of 
the different phases is much smaller than the electronic 
kinetic energy. The most stable phase around a: = is the 
antiferromagnetic phase with orbital ordering. One of the 
three t2g orbitals at the Fe site is occupied, while the two 
others are pushed to higher energies. This phase is insu- 
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FIG. 3: (Color online). Total energy (per site) of the ferri- 
magnetic solutions with and without orbital ordering and for 
the antiferromagnetic solution (orbitally orderered) as a func- 
tion of band filling, x. The splitting between the Fe and Mo 
levels is A = epc — emo + t//3 = — 4fFo_Mo- Inset: band filling 
(vertical axis) as function of chemical potential (horizontal 
axis) in units of tpo-Mo. Note that the doping is equal to 
the band filling minus one. OO stands for orbitally ordered, 
and NO for orbital degeneracy. FM and AF stand for the 
ferrimagnetic and antiferrromagnetic phases discussed in the 
text. 



lating, as shown in the kink of the energy vs. x curve, and 
the finite range of values of the chemical potential com- 
patible with the band filling equal to one (doping a; = 0), 
shown in the inset. This range gives the gap in the elec- 
tronic spectrum, which is of order 1 — 2 x ipo-Mo- Our 
calculations suggest a wide range of phase coexistence, 
between the stoichiometric case, with doping x — 0, and 
doping X ~ 1. At sufficiently high dopings, an orbitally 
ordered ferromagnetic phase is stable. 

The nature of the insulating phase in this regime is 
in agreement with the band structure calculations in2i. 
One of the bands with spin antiparallel to the core Fe 
spin lies separated by a gap from the others. This band 
can accommodate one electron, and doping beyond x ^ 
requires the filling of states above this gap. 

The opposite case, where the Fe and Mo orbitals 
are optimally hybridized (A = 0), is shown in Fig.^. 
This choice of parameters is probably adequate for 
Sr2-a;Laa;FeMo06. The most stable phase near cc = is 
a metallic ferrimagnetic phase, without orbital ordering. 
At sufficiently high dopings, however, the antiferromag- 
netic, orbitally ordered phase prevails. These results are 
consistent with the calculations in^^, although the anti- 
ferromagnetic phase is shifted towards larger values of x, 
due to the different value of iMo-Mo used here. There is 
a first order phase transition between these two phases, 
with a significant range of doping values for which phase 
separation takes place. The two competing phases are 
metallic throughout the physically relevant doping range. 



FIG. 4: (Color online). As in figuregl, but for A = 0. 




FIG. 5: (Color online). As in figure|21, but for A = — Sfpc-Mo. 



The most complex intermediate situation between these 
two regimes takes place for A — — Stpc-Mo, as shown in 
Fig. [3]. The orbitally ordered antiferromagnetic phase 
competes with the ferrimagnetic phase with and without 
orbital order. At low dopings, x < 0, the ferrimagnetic 
phase is the most stable. The insulating, orbitally or- 
dered antiferromagnetic phase prevails near zero doping, 
X = 0, and the ferrimagnetic phase with orbital order has 
the lowest energy for x > 0.6. There are sizable regions 
of phase separation around x — 0. 

The difference between the orbital levels at the Fe ion, 
which can be used to characterize the orbital order is 
shown in Fig.[S]. At low band fillings, when orbital order 
exists, the level splitting is, roughly. Si — 62 ^ Ux. This 
effect leads, for a; = 1 to a sizeable gap between the high- 
est occupied and the lowest empty bands, as discussed 
above. 

The phase diagram which arises from these calculations 
is shown in Fig.0. There is a competition between the 
different phases considered here, with regions of phase 
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FIG. 6: (Color online). Splitting of the orbital levels at the 
Fe ions for A = — 2tFe-Mo as function of band filling, x, in 
units of ipc-Mo. 
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FIG. 7: (Color online). Schematic view of the phase diagram 
as function of the separation between the Fe and Mo(W)levels 
and the band filling x. FM and AF denote the ferrimag- 
netic and antiferromagnetic phases described in the text. OO 
stands for orbital order, and NO labels a phase without orbital 
order. The antiferromagnetic phase at x = 1 is insulating. 
Phase separation takes place in the blank regions. 



structure of Sr2_a:Laa;FeMo06 and Sr2_a:Laa;FeW06. 
The model used resembles the double exchange model 
proposed for the manganites. The presence of a non 
magnetic ion (Mo,W) with levels at similar energies to 
the orbital levels of the magnetic ion (Fe) induces signif- 
icant changes in the resulting magnetic structures. The 
large (threefold) degeneracy of the t2g levels involved also 
influences the phase diagram. 

The presence of non magnetic ion favors phases with- 
out net magnetization, if the band filling is such that the 
conduction electrons tend to reside in the non magnetic 
ion. This situation can be relevant for Sr2-a:Laa;FeMo06 
at sufficiently large dopings. At lower dopings (band 
filling X — > 0), and in a strongly hybridized band, the 
material tends to be a ferromagnetic metal, through a 
mechanism similar to double exchange. 

Orbital order is favored when the conduction elec- 
trons are mostly localized in the magnetic, highly corre- 
lated ion (Fe). Then, a single band splits from the rest, 
leading to a system which is insulating at stoichiome- 
try, a; = (band filling x = 1). This phase can be 
either ferrimagnetic or antiferromagnetic, with a slight 
tendency towards the latter. This situation may describe 
Sr2_a;Laa;FeW06. We have not analyzed the suppres- 
sion of orbital order at finite temperatures. The exis- 
tence of a large gap in the insulating phase adequate for 
Sr2_a;La3;FeW06 implies that, within the Hartree Fock 
approximation, orbital order persists above the Neel tem- 
perature. This insulating phase is separated by regions 
of phase separation from other stable phases at differ- 
ent dopings. A similar situation arises in the single band 
Hubbard model at half filling^. 

The orbital ordering may be difficult to observe in 
a system which tends to be highly disordered, such as 
Sr2_a;Laa;FeMo06 and Sr2_a;Laa;FeW06. Orbital order 
favors anisotropic electronic properties, at each of the 
t2g orbitals selects a plane perpendicular to one of the 
axes of the cubic lattice. The formation of domains be- 
cause of antisites or other types of disorder can restore 
isotropy, although the system should be anisotropic at 
short length scales. 
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